Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics     

A. Mignone  G. Bodo 《Monthly notices of the Royal Astronomical Society》2006,368(3):1040-1054

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.  相似文献   

Angular-diameter distance estimates from the Sunyaev–Zeldovich effect in hydrodynamical cluster simulations     

S. Ameglio  S. Borgani  A. Diaferio  K. Dolag 《Monthly notices of the Royal Astronomical Society》2006,369(3):1459-1468

  相似文献   

Virtual Huanghe River System： Framework and Technology   总被引：2，自引：0，他引：2  

LU Heli LIU Guifang SUN Jiulin 《中国地理科学(英文版)》2006,16(3):255-259

1 Introduction Huanghe (Yellow) River basin is located in 32°–42°N, 96°–119°E. The area of the catchment is more than 752,000km2. The river is 5464km long with a drop in elevation of 4830m. Among the whole area, the moun- tainous and stone area accounts for 29%, loess and hills area 46%, sandy area 11% and plain area 14%, respec- tively. Different natural landscapes exist in this area. The Huanghe River flows through the Loess Plateau, where the soil is eroded seriously (Wang, 2002;…  相似文献   

Radiative Transfer in Inhomogeneous Atmospheres. III     

A. G. Nikoghossian 《Astrophysics》2004,47(3):412-421

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.  相似文献   

应用ATA的对数周期馈源方案于FAST的可能性探讨     

吴盛殷  熊继衮  吕善伟  张江林  徐祥  南仁东  金乘进 《天文学进展》2004,22(2):179-187

介绍应用于美国凤凰计划的Allen望远镜阵(ATA)的对数周期馈源(LPA)概况，初步给出此种馈源的基本参数及性能，仿真估算了其驻波比和方向图，并结合正在进行中的我国500m孔径球面射电望远镜(FAST)计划，探讨了应用这种宽带馈源的可能性及其限制。  相似文献   

A qualitative and quantitative evaluation of experimental model catchment evolution     

G. R. Hancock  G. R. Willgoose 《水文研究》2003,17(12):2347-2363

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

High-resolution calculations of merging neutron stars – II. Neutrino emission     

S. Rosswog  M. Liebendörfer 《Monthly notices of the Royal Astronomical Society》2003,342(3):673-689

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature-dependent equation of state and a multiflavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them with the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos that are produced in the hot, neutron-rich, thick disc of the merger remnant. The cooler parts of this disc contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This imposes constraints on the ability of neutron star mergers to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen 'pole-on' would appear much brighter in neutrinos than a similar event seen 'edge-on'.  相似文献   

The efficiency of the spiral-in of a black hole to the Galactic Centre     

Piero F. Spinnato  Michael Fellhauer  Simon F. Portegies Zwart 《Monthly notices of the Royal Astronomical Society》2003,344(1):22-32

We study the efficiency at which a black hole or dense star cluster spirals in to the Galactic Centre. This process takes place on a dynamical friction time-scale, which depends on the value of the Coulomb logarithm (ln Λ). We determine the accurate value of this parameter using the direct N -body method, a tree algorithm and a particle-mesh technique with up to two million plus one particles. The three different techniques are in excellent agreement. Our measurement for the Coulomb logarithm appears to be independent of the number of particles. We conclude that  ln Λ= 6.6 ± 0.6  for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, such as a dense star cluster, ln Λ is smaller, with a value of the logarithm argument Λ inversely proportional to the object size.  相似文献   

有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献